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ABSTRACT 

The recent analyses of the light curves provided by CoRoT have revealed pul- 
sation spectra of unprecedented richness and precision, in particular, thousands 
of pulsating modes, and a clear distribution of amplitudes with frequency. In 
the community, some scientists have started doubting about the validity of the 
classical tools to analyze these very accurate light curves. 

This work provides the asteroseismic community with answers to this question 
showing that (1) it is physically possible for a star to excite at a time and with 
the observed amplitudes such a large number of modes; and (2) that the kinetic 
energy accumulated in all those modes does not destroy the equilibrium of the 
star. Consequently, mathematical tools presently applied in the analyses of light 
curves can a priori be trusted. This conclusion is even more important now, when 
a large amount of space data coming from Kepler are currently being analyzed. 

The power spectrum of different stellar cases, and the non-adiabatic code 
GraCo have been used to estimate the upper limit of the energy per second 
required to excite all the observed modes, and their total kinetic energy. A 
necessary previous step for this study is to infer the relative radial pulsational 
amplitude from the observed photometric amplitude, scaling our linear pulsa- 
tional solutions to absolute values. The derived upper limits for the required 
pulsational energy were compared with 1) the luminosity of the star; and 2) the 
gravitational energy. We obtained that both upper energy limits are orders of 
magnitude smaller. 

Subject headings: stars: oscillations - stars: individual (HD 174936, HD 49434) - stars: 
variables: 5 Scuti - stars: variables: other 



- 3- 



Introduction 



The first data provided by CoRoT for S Scuti stars made clear tha t the larger th e pre 



cision of the observations, the large r the number of observed modes (IPoretti et al. 



2009 



Garcia-Hernandez et al 



( iGarcia-Hernandez et al 



2009|). In the case, for example, of HD 174936 
20091 ). a 5 Scuti target of the initial CoRoT run, the num- 
ber of frequencies extracted from the observed light curves with the traditional tools can be 
of the order of thousands. The same can be said from other CoRoT 5 Scuti and 7 Doradus 
targets. In such a scenario, the question if all these frequencies are real or not should be a 
previous step to any modeling attempt. To answer this, our aim is to check if a S Scuti or a 
7 Doradus star has enough energy to excite such a large amount of modes simultaneously 
and with the observed amplitudes, for which we need to know the real energy of the modes. 

Up to our knowledge, the only existing studies about the real energy of the modes, 
through characterisation of their amplitudes, are for stochastic pulsators (as solar-like stars) 
which were done through a description of the non-linear effects in the pulsation equations 



fjSamadi et al 



2003 



Belkacem et al. 



20081 ). In the case of thermodynamic pulsators, as 



those driven by a k mechanism, this kind of studies have not yet b een attempt ed, since the 
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1990). Therefore, 



non-adiabatic studies have been done in the linear approximation (lOsakil 
the solutions yielded by the pulsation codes are not absolute, as they are affected by a 
constant indefinite factor, and only ratios of observed and theoretical quantities can be 
compared. 



We sea 



(IMoya et al. 



ed the solutions given by the non-adiabatic linear pulsational code GraCo 



2001 



Moya fc Garrido 



20081 ) to absolute values using the observational 
amplitudes provided by space and ground-base data, as described below. In this way, we 
were able to compare the derived absolute values of the required pulsational energy to 



simultaneously excite all the observed modes, with the real energy available for exciting the 
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modes. 



In our study we use one of the models best repr oducing the ground-based and space 



Garcfa-Hernandez et al. 



2009T). and 



observations obtained for the S Scuti HD 174936 (see 
a model at the centre of the observed HR photometric error box for the hybrid 5 Scuti 
- 7 Doradus star HD 49434. In order to check the consistency of our procedure, we use 
a representative model of the largest amplitude pulsator known, Classical Cepheid stars 
presenting variations of the order of magnitudes, and verified that it is energetically stable 
under its pulsations. 



2. Scaling linear solutions to absolute values 

To determine the energy needed by the system to unstabilize all of the observed modes, 
we have to face the problem that most of the pulsational codes, and in particular the non- 
adiabatic ones, use the linear approximation to solve the system of differential equations. 
Aside from non-linear physical effects, negligible for small pulsational amplitudes, this 
means that every eigenfunction of an oscillation mode is known except for a multiplicative 
constant, i.e., if fi(x) are the solutions of the set of differential equations for a given mode, 
Afi(x), A being a constant, are also solutions. Therefore, before doing any study about the 
global properties of the star, we have to determine the value of this constant (A rea i) scaling 
the linear solutions of each mode to the real values. 

In the linear approximation, the problem of this indefinite constant is solved by scaling 
all the solutions of the pulsation equations to an arbitrary normalization condition. Most 
pulsation codes use as normalization condition 



— = 1 at r = R 

r 



(1) 



- 5 - 



where £ r /r is the radial eigenfunction (IDziembowskil Il97ll ) . The real pulsation amplitude 
on the surface is given by 

And following Eq. [TJ 

5R _ 

— A rea i (3) 

Therefore, the multiplicative constant A rea i, scaling all the linear solutions to the real values 
of the pulsation, may be obtained through the oscillation amplitude of the relative radius 
at the surface of the star. This value can be deduced from the oscillation amplitude of the 
observed luminosity. 



We use the observed amplitude spectrum for HD 174936 (Fig. [Tj and 



Garcfa-Hernandez et al. 



20091 ) as an example to illustrate the procedure here presented. The spectrum shows a 
maximum oscillation amplitude of about 2 mmag, which would correspond to a low-degree 
mode. Decreasing amplitudes are found at both sides of this highest amplitude mode. 
Therefore, we assume that: "The mode amplitude is given by a Gaussian distribution 
centered in the largest amplitude mode, with the maximum amplitude of this mode, and the 
variance fixed by the decrease in amplitude of the modes closest to the highest one. This 
intrinsic amplitude is independent of the spherical degree I of the mode" . The oscillation 
amplitudes of the other modes, out of the highest amplitude one, would be affected by 
geometrical visibility, which is translated in amplitudes lower than the predicted Gaussian 
envelope. 

In this case, the distribution is centered on 377 yuHz, with a maximum amplitude 
of 2.12 mmag and a square root of the variance a = 50 /iHz. The assumed oscillation 
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amplitudes are also shown in Fig. [T] (left, plus signs), compared with the observed amplitude 
spectrum. We notice that a significant number of observed amplitudes fall out of the 
range [200, 500] //Hz, where our distribution predicts amplitudes close to zero. Therefore, 
we built a modified distribution as the addition of the Gaussian already explained plus a 
flat-background at 0.3 mmag for the modes with negligible amplitude (Fig. (TJ right, plus 
signs). There is no physical reasons for this election, but only a fitting of the observations. 
The value of the background was chosen to take into account the amplitude of the observed 
mode at 158.1 yuHz, and the distribution is spanned in the range [100,800] /iHz well 
comprising the canonical range of observed modes for S Scuti. We deliberately leave out 
modes under 100 /iHz as they lay far out of the instability range of 6 Scuti and no model is 
able to predict their excitation. 

The next step is to relate this intrinsic oscillation amplitude in luminosity to oscillation 
amplitude in relative radius. To do so, we use the definition of the apparent magnitude 



variation for a certain wavelength given by 



Balona fc Stobie 



(1197a ) 



Am A = AS X - 2.5 log eAA/A (4) 

where A is the projected area of the photosphere, and AS\ is the variation in surface 
brightness defined as 



AS A = -2.5A(logF A ) (5) 

with F\ being the projected flux at a certain wavelength. 

Given that the unknown multiplicative constant is defined by the variation of the 
radius, we are only interested in the intrinsic amplitude of the radial eigenfunction. 
Therefore, although the relative variation of the projected surface area AA/A includes 
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dependencies as the degree £ of the modes, the viewing angle or the limb darkening, we 
may adopt, for our purposes, AA as the radial variation of the total surface of the star. In 
addition, as the observed CoRoT amplitude spectrum for HD 174936 was taken in white 
light, we will consider F\ as the total flux at the surface. Therefore, Am^ = Am = AM& j, 
since the bolometric magnitude (M& z) and absolute magnitude (m) differ only in a constant. 
Including all this in Eq. 0] yields 



Am = —2.5 loge 



A(lnaT e 4 // )-2— 



(6) 



where a is Stefan Boltzmann constant. Simple calculations give: 

A /? AT 1 

Am = (-5— -10— ^) loge (7) 

Therefore, the variations in brightness are a combination of the variations of the 
relative radius and the relative flux at the surface of the star, measured through relative 
variations of the effective temperature. The observed variations of the relative radius may 
be obtained as 

AR Am , , 



R loge (5 + 10dT) 



where 



dT= 5 Jkl/k at r = R ( 9 ) 
T eff r 

The negative sign of this equation is related to the definition of the stellar magnitude 
and it has no physical origin. We will not take it into account in the rest of this study, 
since it has no influence in the results. Therefore, once the variations in brightness are 



- 8- 



known, the relative radius variations can be obtained, as the variable dT is supplied by the 



non-adiabatic equations. We used the GraCo non-ad iabatic code for t 
ime Dependent Convection description (TDC 



the 



lis purpose, includin g 



Grigahcene et al. 



2005 



Dupret et al. 



20051 ). although its inclusion revealed of no significance in the results. The final relative 
radius variations as a function of the frequency of the modes for the two distributions are 
displayed also in Fig. [TJ 



3. Obtaining the energy balance of the modes for different stellar cases 

Once the multiplicative constant 5R/R scaling all the linear solutions to their real 
values is known, the absolute values of the kinetic energy and the energy required to excite 
the observed modes can be estimated. Then, they can be compared to the energy balancing 
the star and the energy available in the star, respectively. 

To check if the star is able to have thousands of modes excited at a time, we have 
to verify 1) if the energy flux is enough to provide the energy that all these modes need 
per second, and 2) if the total kinetic energy accumulated in them does not destroy the 
equilibrium structure of the star. We note that quantitative estimations of upper limits of 
these energies are enough to achieve the aims of this work. 

Three stellar cases have been studied: 1) a representative model of the classical 
Cepheid stellar type, adopting the largest pulsation amplitude reported (around 2 mag) 
as a falsifiability test of our main assumptions, 2) the 5 Scuti star HD 174936, and 3) the 
hybrid S Scuti - 7 Doradus star HD 49434. 

For the latter two cases, we have used C0R0T and ground-based observations. In the 
case of HD 49434, the publication of the data analysis is still in preparation and we used 
here only some general and preliminary results, which is, however, enough for our purpose. 
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3.1. A classical Cepheid star 



Our first stellar case applies to the variable stellar type with the highest amplitude 
luminosity known to date, the Classical Cepheid type. This type of stars do not show a 
large amount of pulsational modes; however, we are going to use it as a test of falsifiability 
of our main hypotheses, since for them, its high variability offsets the smaller number of 
pulsation modes. 

We built a representa tive equilibrium model lying in the classical Cepheid instability 



strip (ISandage et al 



20091 ). The main physical characteristics of thi s model are 



Bono et al. 



displa yed in 



fll999h . The 



Table 1, and are similar to one of the thoeretical models studied in 
energy balance of the pulsational modes is then compared to some absolute values of the 
star and with some theoretical results obtained through non-linear resolutions, as described 
in Section 2. 

The variation of the relative radius of the fundamental radial mode obtained following 
our procedu re (AR/R = 0.16) is similar to the predicted for an equivalent non-linear model 



described in 



Bono et al. 



fll999h (AR/R = 0.13). 



3.1.1. Pulsational energy 



The energy balance between gains and losses, for a single mode, in a complete period 
of oscillation is given by (jUnno et al.lll989l ) 



W = - / — S(e N - -V • F)olM 



a 



P 



(10) 



This quantity is provided by GraCo and can also be regarded as the required energy 
for each mode to be driven during one oscillation cycle. The scaling of the linear solutions 
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to their real values is done multiplying Eq. 9 by A 2 real . As we are interested in an upper 
limit for this quantity, we assume that for every mode the energy balance is positive, i.e. 
that each mode removes energy from the star. This is only true for unstable modes, as 
stable modes transfer energy to the star. However, in this way, we make sure that the real 
amount of energy subtracted by the modes from the star, regardless of the theory used to 
obtain it, will be always lower than our estimated upper limit. 

The total energy during a cycle of oscillation is given by the sum of all the energy 
interchanged by all the modes in the observed frequency range. We compared this energy 
to the luminosity of the star, as a measure of the available radiation energy per time unit. 
For Classical Cepheid stars, we have used the fundamental radial mode and its overtones. 
Therefore, the sum of W ■ a for each mode (a being the frequency of the mode in Hz) yields 
a total energy per time unit needed by the modes to be overstable of 10 35 ergs s _1 . The 
luminosity of the star is of the order of 10 36 ergs s _1 , that is, one order of magnitude higher. 
As we expected the real energy value to be even lower than the one estimated here, we 
conclude that this extreme case has enough energy per time unit to unstabilize all observed 
modes. This result is used only to legitimize our main assumptions. 



3.1.2. Total kinetic energy 



For the secon d comparison we have calculated the kinetic energy of each mode as 



( jUnno et al. 



1989|) 



l r M - 

Ekin = -O- 2 / \£\ 2 dM 

* Jo 



Once the kinetic energy of each mode is scaled by its variation of the relative radius 
(multiplying by A^ eal ), the total sum of the kinetic energies of the selected modes gives a 
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Bono et al 



(119991 ) give the same order of 



total kinetic energy of the order of 10 41 ergs, 
magnitude for the kinetic energy of the model we are comparing with. This value has to 
be compared also with the total energy of the star, i.e., the sum of the internal energy of 
the star, mass, rotational kinetic energy, turbulent energy, gravitational energy, etc. The 
estimation of all these parameters is a complex task, but the gravitational energy that holds 
the star together is already of the order of 10 49 ergs. Therefore, we can infer that the total 
kinetic energy of pulsations of the radial modes for the case of the Classical Cepheid stars 
is low enough for the star to bear it maintaining the hydrostatic equilibrium. Again, our 
assumptions are validated with this test. 



3.2. A 5 Scuti star observed by CoRoT (HD 174936) 

The 5 Scuti HD 174936 was the first CoRoT target of this stellar type (together with 
HD 50844) presenting unexpected space data. The observational spectrum and the relative 
variation of the radius have already been presented and explained in Section 2 and depicted 
in Fig. 1. The physical characteristics of the model used are displayed in Table 1. 

The total energy per oscillation cycle is given by the sum of all the energy interchanged 
by all the modes in the frequency range [100, 800] /iHz (a generous gamut for 5 Scutis, as we 
said above) and with degrees £ = [0, 7], to be able to account for the thousands of observed 
modes. The observation of s pherical degrees up to I = 7 w i th hig h-precision photometric 



time series were predicted by 



Daszyhska-Daszkiewicz et al. 



(|2006j). 



The sum of W ■ a for each mode yields a total energy per time unit needed by the 
modes to be overstable of 10 29 ergs s , for both proposed amplitude distributions (with and 
without TDC). The luminosity of the star is of the order of 10 34 ergs s _1 , that is, five orders 
of magnitude higher. Consequently, as the real energy value is even lower than the one here 
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estimated, we conclude that the star has enough energy per time unit to unstabilize all the 
observed modes. 

On the other hand, the total sum of the kinetic energies of the selected modes gives a 
total kinetic energy in the range [10 39 , 10 41 ] ergs, depending on the amplitude distribution. 
The gravitational energy holding the star together is of the order of 10 48 ergs. Therefore, as 
in the previous case, the total kinetic energy of pulsations of the thousands of modes is low 
enough for the star to bear it maintaining the hydrostatic equilibrium. 



3.3. An hybrid 7 Doradus - 5 Scuti star observed by C0R0T (HD 49434) 

The last stellar case studied is an hybrid 5 Scuti - 7 Doradus star: HD 49434 which 
has also been observed by C0R0T. The analysis of its light-curve is still in progress and 
the results have not yet been published. Hence, we use only the preliminary result that 
thousands of modes have been detected in the 7 Doradus region, tens of them reported as 
"undoubtly" real (private communication). 

Pulsations in the § S cuti region of the frequency spectrum have also been reported 



( lUytterhoeven et al. 



20081 ) . We assume that in this region the situa tion is similar to tha t 



found for HD 174936 in order to study a possible future situation. 



Uytterhoeven et al. 



( 20081 ) did ground-based multisite observations for this star: we use the highest amplitude 



mode reported, to fix the center and maximum amplitude of our Gaussian distribution. The 
mode has a photometric amplitude of 2.0 mmag at a frequency of 20 /iHz. The observed 
frequencies in the 8 Scuti region present values around 90 /iHz. We have proceed in a way 
similar to the previous case, but using a double Gaussian, one for each region presenting 
pulsational modes. A flat-background of 0.3 mmag has also been used for a second possible 
distribution, as explained in Section 2. 
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We built a model fulfilling the physical characteristics shown in lUytterhoeven et al 



( 120081 ) (see Table 1.). For this study, frequencies in the range [3,900]/iHz, to well cover the 



5 Set and 7 Dor range, and degrees in the range £ = 0,7 have been used. 

We obtain that the total energy per time unit needed by the modes to be overstable 
is of the order of 10 30 ergs s _1 (only using TDC, more suitable for 7 Dor targets). The 
luminosity of the star is of the order of 10 34 ergs s _1 , that is, four orders of magnitude 
higher. We conclude again that the star has enough energy per time unit to unstabilize all 
the possible modes. The total sum of the kinetic energies of the modes selected are in the 
range [10 42 , 10 43 ] ergs, depending on the amplitude distribution. The gravitational energy 
that holds the star together is of the order 10 49 ergs. Again, we conclude that the total 
kinetic energy of pulsation of the thousands of modes is low enough for the star to bear it 
maintaining its hydrostatic equilibrium. 



4. Conclusions 

Using a non-adiabatic linear pulsational code, we have related the observed variation 
in brightness of Classical Cepheid stars, a 5 Scuti, and an hybrid 5 Sct-7 Dor star with 
the intrinsic variation of their relative radius. This has allowed us to scale theoretical 
quantities from this linear code of oscillations to real values of the star, and to make the first 
comparison to date, up to our knowledge, of absolute pulsational quantities with absolute 
parameters of non-stochastic pulsating stars. 

The subsequent analysis of the total kinetic energy of the thousands of observed modes, 
and of the required energy to excite them, shows that all of them have enough energy to 
sustain their oscillation. Our study shows that the results yielded by the analyses of the 
light curves with the classical tools (standard sine wave fitting) dot not come into conflict 
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with what it is physically possible. 
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Fig. 1. — Observed amplitude spectrum of HD 174936 (iGarcia- Hernandez et al.ll2009i ). in- 
trinsic oscillation amplitude (5m, plus signs; scale in the left axis), and relative radius (5R/R, 
crosses; scale in the right axis). Left: built with one distribution. Right: built with two 
distributions (see text for details). 
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Table 1: Physical characteristics of the model used 

Type M/M & L/L Q R/R G 

Class. Ceph. 5.00 496.5 22.49 

5 Set. 1.63 14.21 1.97 

Hybrid 1.60 6.80 1.58 



